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We study the spontaneous symmetry breaking (SSB) of a superfluid Bose-Fermi (BF) mixture in a 
double- well potential (DWP). The mixture is described by the Gross-Pitaevskii equation (GPE) for 
the bosons, coupled to an equation for the order parameter of the Fermi superfluid, which is derived 
from the respective density functional in the unitarity limit (a similar model applies to the BCS 
regime too). Straightforward SSB in the degenerate Fermi gas loaded into a DWP is impossible, as 
it requires an attractive self-interaction, while the intrinsic nonlinearity in the Fermi gas is repulsive. 
Nonetheless, we demonstrate that the symmetry breaking is possible in the mixture with attraction 
between fermions and bosons, like ''"K and *^Rb. Numerical results are represented by dependencies 
of asymmetry parameters for both components on particle numbers of the mixture, A^f and A^b, 
and by phase diagrams in the {Nf, Nb) plane, which displays regions of symmetric and asymmetric 
ground states. The dynamical picture of the SSB, induced by a gradual transformation of the 
single-well potential into the DWP, is reported too. An analytical approximation is proposed for 
the case when GPE for the boson wave function may be treated by means of the Thomas-Fermi (TF) 
approximation. Under a special linear relation between A'f and A^b, the TF approximation allows 
us to reduce the model to a single equation for the fermionic function, which includes competing 
repulsive and attractive nonlinear terms. The latter one directly displays the mechanism of the 
generation of the effective attraction in the Fermi superfluid, mediated by the bosonic component 
of the mixture. 
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I. INTRODUCTION 



The achievement of quantum degeneracy in bosonic 
[l| and fermionic gases of alkali atoms has opened 
the way to the investigation and manipulation of novel 
states of atomic matter, such as Bosc-Einstein conden- 
sates (BEC) d, [3| and superfluid Fermi gases A 
simple but reliable theoretical tool for the study of these 
trapped degenerate gases is the density-functional theory 
0. In particular, the Gross-Pitaevskii equation (GPE), 
which accurately describes BEGs in dilute gases, is the 
Euler-Lagrange equation produced by the Thomas- Fermi 
(TF) density functional which takes into regard the in- 
homogeneity of the condensate In parallel to that, 
many properties of superfluid Fermi gases with balanced 
(equally populated) spin components, and the formation 
of various patterns in them, can be accurately described, 
under the conditions of the BGS-BEG crossover, by an 
extended TF density functional and its time-dependent 
version, as it has been shown recently [rl-fioj. 

One of the fundamental effects in nonlinear media, in- 
cluding BEG; which has been studied in detail, is the 
spontaneous symmetry breaking (SSB) in double- well po- 
tentials (DWPs). Asymmetric states trapped in symmet- 
ric DWPs are generated by symmetry-breaking bifurca- 
tions from obvious symmetric or antisymmetric states, 
in the media with the attractive or repulsive intrinsic 
nonlinearity, respectively [ll[ [the SSB under the action 



of competing attractive (cubic) and repulsive (quintic) 
terms was studied too \XM]- terms of BEG and other 
macroscopic quantum systems, the SSB may be realized 
as a quantum phase transition, which replaces the original 
symmetric ground state by a new asymmetric one, when 
the strength of the self-attractive nonlinearity exceeds a 
certain critical value. Actually, a transition of this type 
was predicted earlier in the classical context, viz., in a 
model of dual-core nonlinear optical fibers with the sclf- 
focusing Kerr nonlinearity [l^. Still earlier, the SSB of 
nonlinear states was studied, in an abstract form, in the 
context of the nonlinear Schrodinger equation (NLSE) 
with a potential term [iJIj as well in the discrete self- 
trapping model [iBl ■ The latter approach to the descrip- 
tion of the SSB effects was later developed in many works 
in the form of the two-mode expansion, each mode rep- 
resenting a mode which is trapped in one of the potential 
wells (see Refs. and references therein). As concerns 
the interpretation of the SSB as the phase transition, it 
may be categorized as belonging to the first or second 
kind (alias sub- or supercritical SSB mode), depending 
on the form of the nonlinearity, spatial dimension, and 
the presence or absence of a periodic external potential 
(an optical lattice) acting along the additional spatial di- 
mension (if any) [l7|. In the experiment, the self-trapping 
of asymmetric states has been demonstrated in the con- 
densate of *^Rb atoms with repulsive interactions [3 • 

Theoretical studies of the SSB in BEGs were extended 
in various directions. In particular, the symmetry break- 
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ing of matter-wave solitons was predicted in various 
two-dimensional (2D) DWP settings [13, including the 
spontaneous breaking of the skew symmetry of solitons 
and solitary vortices trapped in double-layer condensates 
with mutually perpendicular orientations of quasi-one- 
dimensional optical lattices induced in the two layers [l9l | . 
A different variety of the 2D geometry, which gives rise to 
a specific mode of the SSB, is based on a symmetric set 
of four potential wells [1^ (a three-well system was con- 
sidered too [111). Recently, self-trapping of asymmetric 
states was predicted in the model of the BEC of dipo- 
lar atoms, which interact via long-range forces |22|. SSB 
was also studied in the context of the NLSE with a gen- 
eral nonlinearity [l^ . The symmetry breaking is possible 
not only in linear potentials composed of two wells, but 
also in a similarly structured pseudopotential, which is 
produced by a symmetric spatial modulation of the non- 
linearity coefficient, with two sharp maxima f23 |. 

Another generalization is the stud y of the SSB in two- 
[25j and three-component (spinor) |26| BEC mixtures, 
where the asymmetry of the density profiles in the two 
wells is coupled to a difference in distributions of the dif- 
ferent species. Further, the analysis was extended to a 
Bose- Fermi (BF) mixture in Ref. [l^l, where a "frozen" 
fermion component was treated as a source of an addi- 
tional potential for bosons. Dynamical manifestations 
of the symmetry breaking (Josephson oscillations) in a 
Fermi superfluid trapped in the DWP were recently con- 
sidered too [H]. 

In spite of many realizations of the SSB studied in the 
models of degenerate quantum gases, the self-trapping of 
stationary asymmetric states has not yet been considered 
in fermionic systems. An obvious problem is that a Fermi 
gas, loaded into a DWP, cannot feature a direct self- 
attractive nonlinearity, which is necessary to induce the 
SSB in symmetric states. 

The objective of the present work is to introduce a 
model in which the SSB in a trapped Fermi superfluid 
is possible due to an effective attraction mediated by a 
bosonic component, mixed with the fermionic one. Ac- 
tually, we consider the SSB in semi-trapped BF mix- 
tures, with the DWP acting on a single species, either the 
fermionic or bosonic one, as this setting may be sufficient 
to hold the entire mixture in the trapping potential, and 
induce the SSB in its fermionic component. The anal- 
ysis will be performed in the framework of a mean-field 
model, which couples, via nonlinear attraction terms, the 
GPE for the bosonic component to an equation for the 
fermionic order parameter, derived from the respective 
density functional. Inducing an effective boson-mediated 
attraction between the fcrmions requires an attractive BF 
interaction. For this purpose, we take well-known physi- 
cal parameters corresponding to the '^^Rb— ^"^K mixture, 
which features repulsion between rubidium atoms and 
attraction between the rubidium and potassium, charac- 
terized by the respective positive and negative scattering 
lengths, qb > and qbf < [1^. We consider the case 
when the spin-balanced fermionic component of the mix- 



ture is in the unitary regime, corresponding to a diverging 
scattering length which accounts for the interaction be- 
tween the fermionic atoms with opposite orientations of 
the spin, ap — > ±c» (while the BCS regime corresponds 
to the vanishingly weak attraction, with ap — >■ —0; in 
either limit, the effective fermionic Lagrangian does not 
depend explicitly on ap). In fact, the same model with 
a different coefficient of the effective self-repulsion in the 
Fermi superfluid applies to the description of the BF mix- 
ture with the fermionic component falling into the BCS 
regime. Although the self-interaction, induced by the 
quantum pressure, in the equation for the fermionic or- 
der parameter is always repulsive, we demonstrate that 
the SSB in the fermionic component is indeed possi- 
ble in the ^^Rb — K mixture, due to the BF attrac- 
tion which, as said above, mediates an effective attrac- 
tion force in the Fermi superfluid. We also conclude 
that the attraction can induce symmetry-preserving or 
symmetry-breaking localization of both components in 
the semi-trapped mixture, depending on the numbers of 
the bosons and fermions in it. 

The paper is organized as follows. The model is formu- 
lated, in a sufficiently detailed form, in Section II. Results 
produced by the numerical analysis are reported in sec- 
tion III, for two variants of the model, with the DWP 
acting cither only on the fermions, or on the bosons. 
In Section IV, we report approximate analytical results, 
obtained by means of the TF approximation applied to 
the GPE for the bosonic wave function. In particular, 
assuming a specific linear relation between the fermion 
and boson numbers, we can reduce the model to a single 
equation for the fermionic wave function with compet- 
ing self-repulsive and self-attractive terms, the latter one 
explicitly demonstrating the mechanism of the effective 
attraction between the fermions mediated by "enslaved" 
bosons. The analytical results offer a qualitative explana- 
tion to general findings produced by the numerical anal- 
ysis. The paper is concluded by Section V. 



II. THE MODEL 

Our starting point is a model for the degenerate rar- 
efied quantum gas composed of Nb condensed bosons of 
mass jTiB and A^p fermions of mass mp, in two equally- 
populated spin components, at zero temperature. The 
fermionic component is assumed to be in the superfluid 
state at unitarity or, alternatively, in the BCS regime. 
The system is made effectively one-dimensional (ID), as- 
suming that the gas is confined in transverse directions 
by a tight axisymmetric harmonic potential, with trap- 
ping frequencies lo±b, for the bosons and fermions, 
respectively. 

Within the framework of the density-functional theory 
for supcrfluids [sl, Q , the 3D action of the BF mixture is 



{Cb+Cf + Cbf) d-'r dt , 



(1) 
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where £b is the ordinary bosonic Lagrangian density, 



dipB 



^1 



or. 



dt 
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nriB 



(2) 



'0B(r,i) is the macroscopic BEC wave function, and the 
confining potential for the bosons is 



[/B(r) ^ imBwie^^ + ^b(z) , 



(3) 



with R the transverse cyhndric radial coordinate and 
Vb(z) the potential acting in the axial direction 
z. The bosonic superfluid velocity is VB(r,<) = 
(fi/mB)V0B(r, i), where (?B(r,i) is the phase of wave 
function, ■i/'B(r, = ^n-B{v,t)e^^^^^'^\ and nB(r,t) is the 
bosonic density. 

The Galilean-invariant Lagrangian density £f of the 
Ferm i g as with two equally-populated spin components 

is [i,[3q 



£f = i — 



' dt 
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where 'ipp{r,t) is the superfluid order parameter of the 
Fermi gas at unitary [8|, 2mp is the mass of a pair of 
fcrmions with spins up and down, and the potential act- 
ing on the fermionic atoms is 
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UFir) = -rriFLoipR' + VFiz) 
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cf. its bosonic counterpart ([3]). The fermionic su- 
perfluid velocity is VF(r,i) ~ {h/2mp)\79p{r,t), where 
9F{r,t) is the phase of the order parameter, 'ipF{r,t) = 
^nF(r,i)e*^'''^'"'*\ and nF{r,t) is the density of fermionic 
atoms. Constant ^ in expression (|4]) is .J = 1 in the deep 
BCS regime, and ^ ~ 0.4 at the unitarity In calcu- 
lations reported below, we fixed ^ = 0.45, assuming the 
unitarity regime. Lastly, the Lagrangian density £bf in 
expression ([T]) accounting for the BF interaction is 



27rh qbf I , |2 I , ,2 



(6) 



where = mBmF/{mB +mF) is the respective reduced 
mass, and as said above, abf < corresponds to an 
attractive BF interaction, which is necessary to support 
the SSB of the fermionic component in the presence of the 
DWP. We stress our assumptions: i) obf is independent 
of the spin component; ii) the fermionic density profiles 
are identical for the two components. 

In the nearly ID configuration, transverse widths of 
the atomic distributions are determined by the width of 
the ground states of the respective harmonic oscillators: 



a±B 



y^h/{mBUJ±B), a±F = \/h/{2m^ 



(7) 



for the condensed bosons and superfluid fermions at the 
unitarity. The boson and fermion components exhibit 
the effective ID behavior if their chemical potentials are 
much smaller than the corresponding transverse-trapping 
energies, hu!±B and fiwxF/2. Notice that the effective 
mass 2mF in the fermionic harmonic length of Eq. ([7]) 
is a direct consequence of the coefficient 1/8 in gradient 
term jVV'.Fp of the fermionic Lagrangian Q. This coef- 
ficient is chosen to reproduce very accurately the Monte 
Carlo simulations of the 3D unitary Fermi gas [1, [l^] ■ 
Under these conditions, we can adopt the known factor- 
ized ansdtze for the 3D wave functions 13111 . 



V'B(r,<) 
ipF{r,t) 



-J?V(2<b) 



-1/2 



g-i?V(2alF) 

'^F jyjT ^F{z,t) 



a±F 



(8) 
(9) 



where ^b{z) and $f(-z) arc the ID (axial) wave functions, 
which are subject to the usual normalization conditions. 



+ 00 



i{z, i)p dz 



+ 00 



|$F(^,i)l' dz 
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Inserting expressions ([8|) and (HJ into Eq. ([TJ , the action 
can be written as S" S'^^^' - [Nb huj±B - Np huj±p /2] t, 
with the effective ID action 



- / dt 



+ 00 



dz 
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BF 
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where the usual ID Gross-Pitaevskii Lagrangian density 
is 



(ID) 
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dt 
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where the boson self-interaction strength in ID is 

GB^2NBhwi_BaB. (13) 

Further, the ID fermionic Lagrangian density in Eq. pT|) 
is 



,(1D) 



dt 



^ 5*F 



8mp 



dz 
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with the effective strength of the fermionic quantum pres- 



A = (37r2)2/3(3e/5)ft2Ar2/3/(2^p„V3)^ (^5) 

emerging as the coefficient in front of the bulk kinetic 
energy of the Fermi gas in the unitarity limit. Finally, 
the ID Lagrangian density of the BF interaction is 



(ID) 
BF 



-NbNp GbfI^bI'I^fI 



(16) 
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where the corresponding ID interaction strength is 

Gbf = 2fi^aBF/("^Ra_LBa_LF)■ (17) 

For numerical calculations, we set a^B = aj_F = o,± 
and wxB = i^XF = i^x, measuring lengths and time in 
units of ax and w^^, respectively. This implies that 
2mF = mB (hence, mR = tob/3), a condition which 
IS roug hly satisfied by the ^"^Rb Kb mixture. We 
consider ^"K atoms in the two equally-populated hyper- 
fine states \F = 9/2, mp -9/2) and \F = 9/2, m^ = 
— 7/2), and the ^^Rb atoms in the hyperfine state \F = 
2,mp = 2). This mixture is a good candidate for experi- 
mental study of the SSB because the BF scattering length 
is negative, as stressed above, cbf ~ — 284ao, where qq 
is the Bohr radius [l^. Simultaneously, the scattering 
length for collisions between rubidium atoms is positive, 
Ob lOSao [1, [ll,!!^- In principle, the magnetic field 
inducing the Feshbach resonance can affect the values of 
cbf and ob, but we neglect this effect in all the calcula- 
tions reported in the present work. 

The application of the variational procedure to the ef- 
fective action ([TT|) produces a system of coupled NLSEs, 

d 

+.9BFiVF|0Fp] (t>B = ij-(t>B , (18) 



the study of ID gap solitons in the BF mixture in Ref. 
[33I I are also tantamount to Eqs. (fT9|) . up to a difference 
in coefficients. 

To find stationary solutions to Eqs. f ([T5)) and ([T^ . we 
employed a Crank-Nicolson finite-difference scheme for 
simulations of the equations in imaginary time, using the 
Fortran codes provided in Ref. |34| . We employed space 
and time steps Aa; = 0.025 and At = 0.001, and a suffi- 
ciently large number of iterations to ensure the conver- 
gence. The stability of the stationary solutions against 
small perturbations was then tested by simulations in 
real time. 

Due to the attractive character of the BF interac- 
tions (aBF < 0), the true ground state of the 3D mix- 
ture collapses towards energy E = —00. Neverthe- 
less, because of the strong transverse confinement, the 
quasi-lD metastable state has an indefinitely long life- 
time [IHl • Actually, stationary solutions generated by the 
imaginary-time integration represent the ground state of 
the effective one-dimensional BF system based on Eqs. 
(fT8|) and (fT9|) - in the same sense as the famous matter- 
wave solitons realize the ground state of the quasi-lD 
condensate of ^Li atoms ISal. 



III. RESULTS OF THE NUMERICAL ANALYSIS 
A. Axially trapped fermions and free bosons 
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+5BFiVB|0B|^ 

z/aj_, T = LO\t, and 



.gF|0F|^/' + M^F(x) 
i d , 



2 OT 



(19) 



$p, Wb = VB/{fuj^), Wf = Vf/(^j_), and the 
renormalized interaction coefficients are derived from ex- 
pressions (ini), (Hi]), and ([T7)) . 



5b = GB/ifii^±a±) ^2{aB/a±) Nb, 
gBF = GBF/(^i^±a_L) = 6aBF/a_L, 

ffF = A/{aY'nL,^)^{i7rY'm/5)N^^'- (20) 

Note that the rescaled wave functions are subject to the 
same normalization conditions as in Eqs. (jlOp . i.e.. 



+00 



\(t)Bix,T)f dx 



+00 



|0F(x,r)p = 1. (21) 



If condition 2mF — tub docs not hold, the coupled equa- 
tions can be cast in the same form, with a difference 
that coefficient 2mF/mB appears in front of the second 
derivative in Eq. ([T5)) . 

The coupled NLSEs in the form of Eqs. (fT9|) actually 
generalize the static and dynamical equations for BF sys- 
tems which were used in various settings in Refs. (32j . In 
particular, the semi-phenomcnological equations used for 



We start the analysis by considering the configuration 
with the DWP acting solely on the fermionic component: 

Wf {x) = aFX^ -f /3f cxp {~^fx^) , Wb{z)^Q, (22) 

where all constants are positive. An elementary consid- 
eration of this potential demonstrates that it features the 
double- well structure provided that ap < PfIf, with two 
symmetric potential minima located at points 



±7f A/ln(/^F7F/aF)- 



(23) 



We reports results of simulations for ap = 1/2, /3f = 16, 
and 7f — 10, which adequately represents the generic 
situation; in this case, Eq. ([23|) yields Xmin ~ ±0.76. 

The purpose of the analysis is to construct x- 
symmetric and asymmetric ground states of the system, 
varying the control parameters, and identify the respec- 
tive SSB bifurcation, i.e., the transition to the asymmet- 
ric ground state. The asymmetry of its Fermi and Bose 
components is characterized by parameters 

„ _ !^\^FMx)?dx- f_^\(t)F,B{x)\'^dx 

J_oo \4>F,B{x)\^dx 

Recall that the denominator in this expression is actually 
1 for both species, as per Eq. ((2T|) . 

In Fig. [T]we display a set of axial (ID) profiles of the 
densities of both components in the ground state, gener- 
ated by the integration of Eqs. ([T5)) and in imaginary 
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FIG. 1: (Color online). Density profiles of tlie ^'^Rb Kb 
mixture, 4>b{^) a-nd (j>p(x), marked by tlie respective labels, 
in tlie case of the double-well potential (shown in arbitrary 
units by curve "pot") acting on the fermions, while no axial 
potential is applied to the bosons. The three panels differ 
by the number of bosons, A^b, as indicated in each panel. 
Recall that the axial coordinate x is measured in units of 
the transverse-confinement length, a±, while (f>B and (f>F a-re 
displayed in units of aj^. 



time, in the case of the BF attraction and weak repulsion 
between the bosons. The respective values of the inter- 
action coefficients in Eqs. ([T5)) and ([T^ are taken as per 
Eqs. (pn]) . with the above-mentioned values of Obf and 
Ob for the ^^Rb — Kb mixture and a fixed transverse- 
trapping length, a± = 1 /.jm, substituting various val- 
ues of atomic numbers A^b and N-p. The figure clearly 
shows the transition from the symmetric ground state 
to the asymmetric one, with the increase of the number 
of bosons, Nb, i.e., as a matter of fact, the strength of 
both the boson-boson and BF interactions. Accordingly, 
the increase of Nb leads to a stronger overlap between 
the bosons and fermions, and to the SSB, i.e., the transi- 
tion from symmetric ground states to an asymmetric one, 
which happens (simultaneously in both components) be- 
tween TVb = 150 and 170, for a fixed number of fermions, 
Nf = 300. 

Similarly to what is shown in Fig. [TJ the symmetric 
ground state is replaced by an asymmetric one with the 
increase of Np at fixed A'b, as this implies the strength- 
ening of the BF interaction. The summary of the re- 
sults for the SSB in the present setting is provided in 
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FIG. 2: (Color onhne). 

bosons (B) and fermions (F) in the '"Rb — Kb mixture, 
loaded into potential (|22|) . as functions of: (a) the number of 
bosons, Nb , at a fixed number of fermions, Nf = 300; (b) the 
number of fermions, at fixed Ab = 500. 



Fig. [2] by plots of asymmetry parameters ([24| versus 
A^B for fixed Np, and vice versa. Naturally, the SSB 
of both components happens at the same point [for in- 
stance, at A^B ~ 105 in panel (a)]. Nevertheless, the re- 
sulting bosonic asymmetry is essentially stronger [in Fig. 
[IJa) - up to a point, A'b ~ 580, at which both 9b and 9p 
attain values very close to 1, i.e., practically all the atoms 
are collected in a single potential well]. In Fig. [TJb), the 
behavior of the fermionic asymmetry, Qp, is different: it 
jumps to a maximum value at the SSB point, and then 
gradually decreases. This difference between the bosonic 
and fermionic components is natural, as the intrinsic re- 
pulsion in the latter one tends to restore the symmetry 
between the distributions of atoms in the two potential 
wells. 

The results are further summarized in Fig. [31 which 
displays the phase diagram of the mixture. There are 
three regions in the (A'fjA^b) plane: an area where 
the attraction to fermions cannot keep bosonic atoms in 
the trapped state ("free bosons"), the region where the 
bosons and fermions are trapped in the symmetric state, 
with respect to the DWP, and the SSB region, where the 
mixture is trapped in the asymmetric ground state. 

An example of the dynamical development of the SSB 
from an initially symmetric configuration, in the case 
where the ground state is asymmetric, is presented by 
Fig. [3l Initially, the bosons and fermions form a stable 
symmetric bound state, via their mutual attraction, in 
the single-well potential acting on the fermions, which 
is taken in the form of dm with Pp = 0. Then, /3f 
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FIG. 3: (Color online). The phase diagram of the *''Rb-'*''Kb FIG. 5: (Color online). The phase diagram of the ^''Rb -° K 
mixture in potential (|22p . which acts only on the fermions. mixture loaded into potential (|25|) which acts only on the 
The diagram shows, on the logarithmic scales in the (A''b, A^f) bosons, 
plane, regions of the symmetric ground state, and of the spon- 
taneous symmetry breaking (SSB). "Free bosons" implies de- 
localization of the bosonic wave function. 

B. Axially trapped bosons and free fermions 





Now, we consider the action of the DWP on the bosons 
only, taking the potential as 
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FIG. 4: (Color online). Contour plots for the evolution of 
densities |(^B(a::,t)P and \<f)-p{x,t)\'^ of the (a) bosons and (b) 
fermions in the *^Rb — Kb mixture. At f = 0, the fermions 
are trapped in the single-well potential (|22|) . with /3f = 0. 
Then, /3f linearly increases from to 16 by t = 80, which 
implies the transition to the double-well potential with well- 
separated symmetric minima. Numbers of particles are A^b = 
500 and A^f = 300. 



is ramped linearly in time (0 < t < 80) from /3f = to 
/3f = 16, which leads to splitting the single potential well 
into two, as per Eq. ([^5]) . The dynamical picture clearly 
shows the transition of the initial symmetric state into 
the symmetry-broken one. Both components get spon- 
taneously collected in one of the wells, where they stay 
together due to the mutual attraction, approaching an 
equilibrium configuration. 



Wb{x) = asx^ + Pb exp (-7Ba;') , Wp = (25) 



For this setting, numerical results are presented with 
Q!B = 1/2, /3b = 16, and 7b = 10, i.e., the same pa- 
rameters of the DWP as used above for the trapping the 
fermions. We again aim to construct the ground state of 
the system as a function of the control parameters, and 
investigate its spontaneous transition into an asymmetric 
shape. 

Results of the analysis for this setting are summarized 
in the respective phase diagram of the mixture plotted in 
Fig. [5] (cf. Fig. In this case too, there are three re- 
gions in the {Np,Nb) plane: an area where the fermions 
cannot be held in a localized state by the attraction to 
bosons ("free fermions"), the region where the fermions 
are trapped, along with the bosons, in a symmetric 
ground state, and the region where the trapped ground 
state is asymmetric ("SSB"), for both the fermionic and 
bosonic components. 

An example of the transition from the symmetric 
ground state of the BF mixture to an asymmetric one, 
caused by the increase of the number of fermions from 
A^F = 10 to A^F = 1000, while the number of bosons 
trapped in potential (|25p is kept constant, is displayed 
in Fig. [6] Although the applicability of the functional- 
density description for A^f = 10 [in panel (a)] may be 
disputed, this figure adequately shows the transition to 
the SSB. 
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FIG. 6: (Color online). An example of the spontaneous sym- 
metry breaking in the ^^Rb-*'^K mixture trapped in potential 
(|25|l . The two panels differ by the number of fermions: (a) 
iVp = 10; (b) Nf = 1000. 



The case of the fermionic component in the 
BCS regime 



The analysis presented above pertained to the super- 
fluid BF mixture with the spin-balanced fermion com- 
ponents in the unitarity regime, where the s-wave scat- 
tering length, ap, which accounts for the interaction be- 
tween the spin-up and spin-down fermions is extremely 
large (ideally, ap ±oo). Actually, essentially the same 
Lagrangian (|4]) also applies to the BF mixture with the 
fermionic component kept in the BCS regime, with ap is 
negative and small (ideally, ap — > —0). In this regime, 
one is practically dealing with a gas of ideal fermions, 
because the respective superfluid energy gap is exponen- 
tially small. In practical terms, to study the system 
whose Fermi component falls into the BCS regime, it 
is sufficient to set ^ = 1 instead of ^ = 0.45 in Eq. ^ 
. Obviously, this change of ^ will make the fermions 
effectively more repulsive [see Eq. p^ ]. as the Pauli re- 
pulsion attains its maximum at ^ = 1. We have verified 
that, as the SSB in the mixture emerges chiefly due to 
BF attraction, the increase of the intrinsic repulsion in 
the Fermionic component corresponding to ^ = 1 makes 
the natural BF attraction in the ®^Rb — '"^ K mixture 
insufficient for the appearance of the SSB for relatively 
small values of iVp and Nb considered above. To achieve 
the transition to asymmetric ground states (for the same 
value of = 1 //m as taken above), it is necessary to 
consider values of A^b and Np exceeding 1500 [cf. Figs. [3] 
and O which display the phase diagrams of the mixture 
with the fermionic component in the unitarity regime for 
A^B,A^F < 1000]. 



IV. THE ANALYTICAL APPROACH: THE 
THOMAS-FERMI APPROXIMATION FOR THE 
BOSE COMPONENT 

A. The general case 

For the application of the analytical approach, we 
use the stationary version of general equations ^TE\\ 
and (Uni); obtained by the substitution of (f>B.F{x,T) ~ 
cxp (— i/iB.F''') UB.p(a;), where chemical potentials for the 
localized states must be non-positive, yUe.F < 0, and real 
functions ub,f obey the following equations (with the 
prime standing for d/dx): 



1 



1 



[-A^B + VFb(.t) + 53^1 + ^BF^^FWp] Ub 



= 0, 
(26) 



- -u'' 
8 ^ 



-^p 



4/3 
5Fltp 



+ 5BF^BU| + Wp(x) 



Up = 0. 
(27) 

An essential simplification of Eqs. ([26|) and (|27l) can be 
achieved if the TF approximation may be applied to the 
former equation, i.e., the term with the second derivative 
may be neglected in it 0, 0, In the present setting, 
with the characteristic size of the trapped states Ax ^ 1 
[see Eqs. ([231) and Figs. [T]and[n], and the wave functions 
subject to normalization conditions (|2ip [hence the am- 
plitudes of normalized density Ug(a;) and w|(x) are also 
~ 1] , a straightforward consideration of Eq. ((26|) demon- 
strates that the kinetic energy (the second derivative) is 
negligible in comparison with either nonlinear term un- 
der conditions TVb a_L/aB or iVp ^ a^/ (10 |aBp|). For 
the value of = 1 //m adopted above, and the values 
of a-Q and jaBpl for the ®^Rb — K mixture, these con- 
ditions reduce to quite realistic inequalities, N-q ^ 200, 
iVp > 10. 

If the TF approximation is valid, it allows one to solve 
Eq. (|26p in the following form: 



„2('^^ _ / [\9'bf\NpuI{x) - WYi{x) mb] , at \x\ < xq, 
^^"^^ ~{ 0, at |x| > xo, 

(28) 

where it is taken into account that we are dealing with 
5bf < 0, and xq is a positive root of equation 



uUxo) ^ i\gBF\ Nf) M^^sl^^o) - Mb] 



(29) 



Note that the bosonic chemical potential, ^Bi is not an 
arbitrary parameter; instead, it must be found from the 
normalization condition (j2ip . applied to expression 



2 / ° [\9bf\ NfuUx) - Wsix) + Mb] dx = ge- (30) 
Jo 



B. A tractable example 

The substitution of approximation 
equation (j27p for the fermionic function 



for Uq into 
up(x), allows 
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one to reduce the underlying system to a single equa- 
tion for u-p(x); however, in the general case this equation 
is quite complex. In particular, the additional equation 
(|29p for xq actually makes the resulting equation for wp 
nonlocal. Thus, in the general case the TF approximation 
does not yield an explicit analytical solution. Neverthe- 
less, it can be obtained in a special case, when Wb = 
[cf. Eq. dH])] and fiB = 0. In this case, Eq. ^ yields a 
simple local relation, 

Ub{x) = (|.gBF| Nf/ob) ul{x), (31) 

which is valid at all x, and Eq. (pO| reduces to a special 
relation between the boson and fcrmion numbers, 

= 9b/ \9bf\ = - (flB/SaBp) Nb, (32) 

where we have made use of Eqs. (PU)) : note that this 
approximation is meaningful only in the case when the 
signs of Ob and obf are opposite. For the parameters 
of the ^^Rb -"^ K mixture Eq. ([32]) amounts to Np ~ 
0.127iVB . Finally, the single equation for the fermionic 
stationary function takes the following form, upon the 
substitution of expression (j3ip : 

-iltp + [iyF(2;) - Mf] Up 
o 

, /q 2\2/3 3^ 2/3 7/3 l^alp Ar 3 n 

+ l''^ J ^^F "f —NpUp = 0, (33) 

5 aBO,± 

where we have again used Eqs. (|20p . Note that the last 
term in Eq. ([33|) directly illustrates the possibility pro- 
posed in this work, namely, that the interaction mediated 
by the boson field may give rise to an effective attraction 
in the Fermi component of the BE superfluid mixture: 
indeed, the coefhcient in front of this term is propor- 
tional to the square of the scattering length, a|p, which 
accounts for the BE attraction. It is also worthy to note 
that, in the present simplest approximation, which makes 
it possible to eliminate the boson field and reduce the 
model to the single equation for the fermionic function, 
the attractive character of the resulting boson-mediated 
interaction requires ob > 0, i.e., the repulsive character 
of the direct interaction between the bosons. Note that 
in the case when both the BE and boson-boson interac- 
tions are attractive, i.e., both ob and gbf are negative, 
the present approximation is impossible, according to Eq. 
([5^ : in the case of the BE repulsion and boson-boson at- 
traction, i.e., Ob < and cbf > 0, the approximation 
is possible, but it leads to an effective boson-mediated 
repulsion between the fermions. 

Equation p3|) is a variant of the stationary NLSE with 
two competing nonlinear terms, the self-repulsive one ^ 
Np^^Up^^, and the self-attractive cubic term. For a± ~ 
1 /um and the scattering lengths corresponding to the 
^^Rb K mixture, the coefficient in front of the cubic 
term is 18a|p/ (obOj,) ~ 0.71, while, with ^ = 0.4 (recaU 
it corresponds to the unitarity regime for the fermion 
component), the coefficient in front of the repulsive term 
is (37r2)2/3(3^/5)si 2.30. 



The SSB controlled by competing nonlinearities, viz.^ 
self-focusing cubic and self-defocusing quintic terms, was 
studied in Rcfs. [l2|, where it was concluded that the 
respective SSB diagrams, showing the asymmetry versus 
the total norm of the mode (cf. Fig. [2]), tend to form a 
closed loop connecting an initial symmetry-breaking bi- 
furcation and a final symmetry-restoring one. The differ- 
ence of Eq. ([55)1 is that here the self- focusing (cubic) term 
has a higher nonlinearity power than its self-defocusing 
counterpart of power 7/3, therefore no closed-loop bifur- 
cation diagram is expected in the present case. 

The substitution of uf{x) = vp/y/Np casts Eq. ((33)) 
in a parameter-free form [the respective equation for vf 
seems as Eq. f (|55)) with Np replaced by 1]; of course, this 
substitution changes normalization (j2ip for the fermionic 
function, as the norm of vp is exactly N-p, but not 1. 
Thus, in the present approximation, the BE mixture is 
described by the universal equation. This circumstance, 
along with condition (j32]) necessary for the applicability 
of Eq. p3p , may explain the fact that the border of the 
trapped states in Fig. [3]is practically a straight line with 
the slope equal to 1, which runs into the SSB border at a 
critical value of Np [the latter one actually corresponds 
to the critical norm of field vp (x) at which the SSB occurs 
in the framework of Eq. ([55)) ]. 

The actual SSB point generated by Eq. ()33p can be 
predicted in an approximate analytical form by means 
of the two-mode expansion, which, as mentioned above, 
is commonly used for the description of SSB effects in 
DWP settings [l^ , assuming that the stationary solution 
is approximated by a superposition of two linear modes, 
u±{x), which are trapped in the left and right potential 
wells: 

uf{x) = A+u+{x) + A_u-{x). (34) 

Symmetric states correspond to = A- = Aq in Eq. 
(|M)) . The SSB border can be found by looking for a point 
where a solution with an infinitely small antisymmetric 
perturbation (6), A± = Aq zt 6, branches off from the 
parent symmetric state. By performing this analysis (we 
do not display straightforward details here), one finds 
that the value of Aq for SSB scales inversely proportional 
to \/Nf, hence the increase of the number of particles 
facilitates the transition to the asymmetric ground state, 
as one may expect. 



V. CONCLUSION 

The objective of this work is to extend the analy- 
sis of the SSB (spontaneously symmetry breaking) in 
DWP settings (double- well potentials), which was re- 
cently studied in BEC and bosonic mixtures, to the BE 
(Bose-Fermi) mixtures. The system is described by the 
GPE (Gross-Pitaevskii equation) for the bosons, which 
is nonlinearly coupled to the equation for the fermionic 
order parameter derived from the density functional in 
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the unitarity limit (in fact, a similar model also ap- 
plies to the BF mixture with the fermionic component 
kept in the BCS regime). Direct symmetry breaking 
in the Fermi superfluid trapped in the DWP is impos- 
sible, as it must be induced by attractive interactions, 
while density perturbations in the degenerate fermionic 
gas interact repulsively. Nevertheless, we have demon- 
strated that the SSB is possible in the mixture of ^^Rb 
and atoms, due to the attraction between fermions 
and bosons. The most interesting case, when the effec- 
tive SSB in the fermionic component could be studied in 
the "pure form" , is that when the fermions are subject 
to the action of the DWP, while there is no potential 
confining the bosons. We have also investigated the al- 
ternative situation, with the DWP acting solely on the 
bosons. Our phase diagrams in the (iVp, A^b) plane, pro- 
duced by means of numerical methods, clearly show that 
the inter-atomic attractions can produce both symmet- 
ric and symmetry-broken localization of the atoms which 
are not subject to the direct action of the trapping po- 
tential. Applying the TF (Thomas-Fermi) approxima- 
tion to the bosonic equation, we have also developed an 



analytical approximation, which allows us to reduce the 
model to the single equation for the fermionic function. 
In the latter case, the model explicitly demonstrates the 
generation of the effective attraction between fermions 
mediated by the bosons. 

The analysis reported in this work can be extended in 
other directions. A straightforward generalization may 
deal with the system including the confining potential in 
both components, as well as a more general analysis of 
the TF approximation. A challenging possibility is to 
predict similar effects in multi-dimensional Bose-Fermi 
mixtures. 
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